Magnetic-Field Generation by Randomly Forced Shearing Waves 
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A rigorous theory for the generation of a large-scale magnetic field by random nonhelically forced 
motions of a conducting fluid combined with a linear shear is presented in the analytically tractable 
limit of Rm <^ Re <g 1. This is a minimal proof-of-concept calculation aiming to put the shear 
dynamo, a new effect recently reported in a number of numerical experiments, on a firm physical 
and analytical footing. Numerically observed scalings of the wavenumber and growth rate of the 
fastest growing mode, previously not understood, are derived analytically. The simplicity of the 
model suggests that shear dynamo may be a generic property of shear flows — with ubiquitous 
relevance to astrophysical systems. 

PACS numbers: 47.27. W-, 47.65.Md, 52.30.Cv, 95.30.Qd 



CI,. 

I ■ 

o 



> 
in 



OO 

o 



X 



Introduction. The magnetogenesis, or origin of cos- 
mic magnetism, is one of the fundamental problems in 
theoretical astrophysics. It has long been believed that 
the magnetic fields observed in most astrophysical bodies 
owe their existence to the dynamo effect associated with 
the turbulence of the constituent plasmas. It is not con- 
troversial that turbulence of a conducting fluid amplifies 
magnetic fluctuations at scales comparable to or smaller 
than the scale of the motions. Small-scale magnetic fluc- 
tuations are indeed observed ubiquitously, but in most 
astrophysical systems, one also finds magnetic fields co- 
herent on scales substantially larger than the scale of 
the turbulence (e.g., Q). Generation of such fields, or 
mean-field dynamo, is expected to require a combined 
action of turbulence and some large-scale-coherent fea- 
ture. One well-established such additional ingredient is 
net helicity (or, more generally, mirror-non-invariance) of 
the motion. The presence of helicity engenders growth of 
large-scale ("mean") magnetic field, known as the a ef- 
fect. While deriving the a effect analytically for realistic 
turbulent systems requires rather drastic closure assump- 
tions, which usually cannot be justified rigorously, it is 
at least certain that the effect exists in certain physically 
meaningful and rigorously treatable limits: for example, 
in chaotic flows of a conducting fluid at low Rm 0, Q . 
The combination of such proof-of-concept analytical re- 
sults, intuitive explanations of how the effect works Q 
and a body of numerical evidence [^, 0] has helped build 
a case for the a effect as a real physical phenomenon. 

It has been suggested 0, H, that even in the absence 
of helicity, mean-field dynamo is possible if a large-scale 
velocity shear is present. The importance of such a pos- 
sibility can hardly be overestimated, as shear is a ubiqui- 
tous feature in astrophysics (usually associated with dif- 
ferential rotation). A recent numerical study showed 



that shear dynamo does exist. However, its nature has re- 
mained poorly understood. The uncertainty is increased 
by the fact that, while the original derivation of the effect 
relied on a quantitative outcome of a closure calculation 
01, the effect proved difficult to identify unambiguously 
by numerical computation of the mcan-field-theory coeffi- 
cients and appeared to go away in rigorously solvable 
limits: the white-noise- velocity model and low-Rm mean- 
field magnctohydrodynamics [1, [l^, (although there 
has been some disagreement on whether shear dynamo is 
possible for Rm < Re < 1 0, Q). 

In this Letter our aim is a minimal proof-of-concept 
calculation that would put the shear dynamo effect on a 
firm physical footing akin to that enjoyed by the a effect, 
we propose a very simple mean-field theory that rigor- 
ously predicts a large-scale dynamo driven by randomly 
forced shearing waves in the limit of Rm ^ Re ^ 1 . The 
effect does not require any adjustable parameters. We 
also recover the scalings of the characteristic wavenum- 
ber and growth rate of the fastest-growing mode that 
were observed in a number of numerical studies 13, 15 1 
but have not so far been explained analytically. 

Shearing Waves. First let us introduce a model veloc- 
ity field that will be used to obtain a dynamo. Consider 
an incompressible fiuid with an imposed background lin- 
ear shear, U = Sxy, and assume that the magnetic field 
is dynamically weak, so the Lorcntz force is negligible. 
Then the velocity deviation from U satisfies 



dtu^ 



SxdyU - 



Su^y + u ■ Vtt = — Vp- 



/, (1) 



where p is pressure determined from incompressibility V • 
u = 0, V is viscosity and / is a random body force, which 
is assumed to have a characteristic scale I / . 

Let us now make two simplifying assumptions. First, 
let Re ~ ulf/v ^ 1, so we can neglect the nonlinear 
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term in Eq. ([T]). Second, let dzU = and dzf = 0, i.e., 
consider a "quasi-2D" velocity, which has all three vec- 
tor components but no z dependence. This velocity will 
help make analytical calculations particularly transpar- 
ent, while, as indicated by numerical experiments [16[, it 
is a favorable but not a uniquely special case vis-a-vis its 
dynamo properties. Since now V • = 0, the xy-plane 
velocity has a stream function: = zx V<i>. Similarly, 
/-'- = z X VF. We seek the solutions of the (linearized) 
Eq. Il]) as superpositions of "shearing waves" : 

^ = J2^'^oit)e"'^'^\ ^..=^^..A,o(^)e"=(*)•^ (2) 

ka feo 

where fco = {kxO,ky,0), k{t) = {k^o — Skyt,ky,0). Then 
the amplitudes of the shearing waves satisfy 

dt [k\t)<i>ko] = -J^fc^(t)$fco + k^it)Fk„, (3) 

dtUzko = -'^k'^{t)Uzka + fzko- (4) 

Eq. ^ was obtained by taking i • [V x Eq. Eqs. 
^ and (dl have explicit solutions 

1>fcoW= f dt'GAt,t')^Fkoin (5) 



p(i) 



Uzko{t)= / dt'G,it,t')fzkoit'), 
Jo 



(6) 



where Gu{t,t') = cxp 



■!l,dt"k\t")\. 

Random Forcing and the Velocity Correlators. In the 
above, we assumed that the body force can be decom- 
posed into shearing waves similarly to the velocity field, 
a convenient modeling choice and a standard approach in 
spectral shearing-box simulations [lo| . However, physi- 
cally, we want to model energy injection occurring at 
some specified characteristic scales (wavenumbers) in 
the laboratory frame, not in the shearing frame where 
wavenumbers depend on time. To implement this, we 
choose the body force to be a Gaussian random field with 

{Fko{t)F*^o{t')) - [AJk'it)] Ak^mt - t'), (7) 
(/.feoW/:feo(^')) = AzAk^mt - t'), (8) 

where Akg{t) = H^k^f - |fcx(i)|)<5|fe„|,fc„j fixes the forcing 
wavenumber kyf = 2iT/lf in the y direction and, via the 
the Heaviside function H, sets the maximum wavenum- 
ber kxf in the x direction. This model is chosen for its 
simplicity; in general, Afc^ (t) can be arranged to give the 
forcing any desired spectrum (e.g., a spherical shell in the 
wavenumber space [lO|). 

A white-in-time in time forcing is a popular choice 
in turbulence modeling 0, [l^l because the amplitude 
of such a forcing can be directly related to the mean 
power injected into the system per unit volume: denot- 
ing e = {u- f) and assuming for simplicity that the same 
power is injected into each velocity component, it is easy 



to show that A± = 2Az = {2Tr/Lxkxf)s/3, where is 
the system size. Note that the white-noise forcing model 
is not particularly restrictive because the resulting veloc- 
ity [Eqs. (O and ^] does have a finite correlation time Tc, 
in our case determined by viscosity, Tc ^ I'j/v (in a more 
general case when Re is not small, the velocity correla- 
tion time is set by the nonlinear terms, so Tc If/urms 
is the typical turnover time of the turbulence; nonrig- 
orously, this case is included in our consideration if the 
viscosity is replaced by turbulent viscosity, v ~ Mrms^/)- 
We shall see that, unlike in the case of other solvable dy- 
namos 0, J it will be essential for the shear dynamo 
that the velocity correlation time is not zero. 
Using Eqs. (I5][5]), the velocity correlators are 

{\<fkoit)\') = A^ fdt'[GAt,t')]'^Akait'), (9) 



(10) 



{\uzkoit)f) ^ Az / dt'[GAt,t')YAko{t'), 
Jo 

while for two-point correlators with t > t' , we have 



{^komi',{t'))^5koKGAt,t')^{\'^ko{t')?). 

(11) 

(".fco(0<fe^(<')) =4„,fe;G.(<,f')(|w.fe„(i')l')- (12) 

Note that the mean hclicity of our velocity field is 
(m • [V X u]) = -2J2ko '^^(*) ("^fco^fco)- ^^"^ ensure 
that it is zero by stipulating {fzk„{t)F^ (t')) = 0. This 
removes the possibility of the a effect. 

Mean-Field Theory. The evolution equation for the 
magnetic field B in the presence of linear shear is 

dtB + SxdyB + u-VB = B-Vu + SBxV + iiV^B, (13) 

where 77 is the magnetic diffusivity. Since the veloc- 
ity field is idependent of z, we can separate the depen- 
dence of B on z by expanding it into Fourier modes, 
B = B{kz)e'^^'^ . Furthermore, only the projec- 
tion B^ onto xy-plane needs to be calculated because 
Bz = {i/kz)'V ■ B-^. For each kz, B^ will satisfy a closed 
equation with kz entering as a parameter and no mode 
coupling in kz. We now seek the solutions of this equa- 
tion again in the form of superposition of shearing waves, 
B-^ = '^ka -^fco (^)6''°^*'' where the perpendicular wave 
numbers fco and kit) are defined in the same way as in 

fco 



the velocity decomposition [Eq. Bu. satisfies 



dtB^^ = SBxkoV - V [k\t) + kl] Bl 

+ Y.^k'„[z^k'{t)]-[k{t)l^ik'{t) 
fcp 

- ik^y^^Uzk'„B^^_^,^, 



B 



ko-k'g 

(14) 
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where I is a unit dyadic. We can separate the large-scale 
field by considering the ko = component of Eq. ((T4)) : 



dtBo = SB.My - V^IBq - ik^ ^ u^koB^ 

ko 



(15) 



(note that B^o = because V • Bq ~ ikzE^o = 0). 

We now calculate B:^^.^^ in Eq. p5|) in terms of Bq via 
Eq. (|14p . This is particularly easy in the limit Rm ^ 
Re <C 1 and S'(///itrms)Rm ^ 1, where Rm ~ u^mslf/v- 
Then 77 A:^ > fc^*, z^fc^^ 5 and |So| » l-B^J, so the dom- 
inant terms in Eq. are r}k^{t)B^ and the k' = k 
components of the wavenumber sums. Then 



B 



-fen 



»^.<feo-Bo + 'i>fcJixfc(t)]fc(t).Bo 



(16) 



r7[/c2(0 + fc2] 

Substituting this into Eq. (fT5|) . we get 

atSo S'Bxoy-[T/ + /3(t)]A:fBo-l-t/c.-ixQ:(t).Bo, (17) 

where /3(0 = |«,fc„(OI V'?[fc'(i) + fc'] « ^7 (negligible 
"turbulent diffusivity" in the limit of low Rm) and 

= E fI4^^S ° MW). (18) 

This can be viewed as a tensor generalization of a familiar 
a term but Eq. (ini) is still a stochastic equation: a(i) 
fluctuates with the correlation time Tc of the velocity field 
and (a) = (cf. [ll|, Q [H ) . 

Averaging Eq. (|17p (with respect to forcing realiza- 
tions), assuming additionally that Stc ^ 1, r/k'^Tc ^ 1, 
and using the standard iterative expansion [20| to calcu- 
late {a{t) ■ Bo{t)), we find the equation for B — (Bq): 



dtB 



Syx - T]kll -klD 



B, 



(19) 



where the term originating from a now looks more like 
a tensor diffusivity (although not sign-definite): 



D 



dt'{[z X a{t)] ■ [z X a{t')]), 



(20) 



where St 1, but t ^ Tc. Substituting for a from Eq. 
(fTS)) . using the fact that and <!> are Gaussian and un- 
correlated with each other (so there is no mean helicity), 
and employing Eqs. pTjl and p^ . we find 



(|^.feo(^OP)(|<&feo(tOP) 
772 [/fc2(t) + fc2] [p(i/)4.;,2] 



(21) 



where K(t, <') = [z x k{t)]k{t) ■ [z x k(t')]k(t'). With- 
out shear, fc(i) = k{t') and K(t,t') = 0. In the presence 
of shear, the wave vectors evolve and we get 



K{t, t') ^{t- t')Skl 



kykx{t ) 

x{t)kx{t ) kykxit) 



-k^ 



(22) 



All matrix elements of D are small compared with other 
coefficients in Eq. (fTO]) and only D^y need be kept because 
there are no other contributions to the xy matrix element 
in Eq. ([1^)) . Eq. now takes the following form 



dtB^ 



-rjkl 
S 



~r]k: 



B 



(23) 



Since K{t, t) — 0, it is manifest in Eqs. (PT|) and 
that D ^ only if the velocity field has a non-zero corre- 
lation time (cf. (rj-fl^). Since we are working in the limit 
Stc <^ 1, the time correlation function of the velocity is 
G'i.(t, t') « cxp [-vk^{t){t - t')] and we have 



Dxy ~ S ^ ky 



ko 



(Kkoimn'fkoitr) 

r]^i^^k»{t) 



(24) 



Since DxyS < is unambiguously sign-definite, there is 
a positive feedback in the dynamo loop in Eq. (p3)) . The 
growth rate is 7 = kz\DxyS\'^/'^ — r/k^. The wavenumber 
and the growth rate of the fastest-growing mode are 



ipeak 



2V 



In 



\D.vS\ 



(25) 



Dependence on Shear. How k^^^^ and 7max depend 
on S* is a key (numerically) measurable parameter de- 
pendence and we should be able to predict it if we are to 
claim that our model calculation contains some universal 
features. The evidence so far has contradicted theory: 
mean-field theories 0, Q predict k^'^'^^ cx S, 7max oc S'2, 



[I0|,ll5| 



while numerical simulations in various contexts 
report kf """^ cx S^/^, 7max oc S. 

Let us substitute into Eq. the velocity spectra 

[Eqs. ([9|) and pO|) ] of the shearing waves excited by our 
model forcing [Eqs. ^ and ([8])]. We get 



(26) 



where 



/ = — 



+ 00 



d^ 



rdre-2-(i+«''/^=H(e/-|e + ^T|) 
Jo 



(27) 



To obtain this formula, we summed over ky, used the 



fact that the integrand in Eq. p4|) is even in ko, replaced 
= {Lx/2ti) J dkxo, changed the integration vari- 
ables to ^ = {kxo — Skyft)/kyf and t = t — t' , denoted 



4 



if — kxf/kyf, Tc 



and Pm = v/r] <C 1. We 



assume S' > 0; if S* < 0, change ky 
wavenumber integral in Eq. 



-ky inside in the 



It is not hard to see that D^y is largest for small 
values of fc^,/, i.e., shearing waves forced in a narrow 
band in kx{t) are particularly effective in amplifying 
the magnetic field [22]. So let us take k^f ~ 2t: / Lx 
(the smallest possible) and ^/ = k^f/kyf <^ 1. In 
this limit the integral (|?7|) is easy to calculate: / = 
(t2/2) [1 - (S'tc/4^/) (1 - e-^^i/^''-)] . This formula en- 
codes the fact that D^y depends crucially on the widths 
of the velocity correlation function and of the spectrum 
of the shearing waves: the Heaviside function in Eq. ([27|l . 
taken as a function of r, has the width / S, which is to 
be compared with Tc- 

If S'tc < < 1 ("weak shear"), we have / « 
whence D^y ~ SPm^uf^^r^kyj/lS, where u^„,s = stc 
(this is clear dimensionally and also easy to show quan- 
titatively for ^ 1). Substituting into Eq. (pS)) gives 



^pcak _ g 



Rm^ 
6^2 



Rm'^Rc 
72 ' 



(28) 



where Re = Urais/vkyf, Rm = Ur^^^/vkyf. These are the 
previously predicted scalings [3,[3|- 

In the opposite limit, ^/ ^ Stc <C 1 ( "not-so-weak 
shear"), we have / w rdf/S. This gives D^y ~ 
Pm'^uf^gT^kyf{2T:/Lx)/9 and, therefore. 



ipoak 




36 



Lxkyf 

(29) 

Such scahngs, k^''^^ cx 5^/^ and 7max « S, have 
been reported in several recent numerical experiments 
lol . . While it is encouraging that they can be recov- 
ered analytically, one should keep in mind that most of 
these simulations were not done in the asymptotic regime 
Rm <C Re <C 1 or at particularly small Stc- The import 
of our derivation is in demonstrating that the scalings 
depend on how strong the shear is and whether the ve- 
locity field decorrelates faster or slower than it takes the 
shearing waves to transit the wavenumber band (in the 
laboratory frame) where they are forced. 

Note that the system size enters explicitly in Eq. 
(PS)) and in Eq. (^5]) [in Eq. , it disappeared because 
we maximized the dynamo by taking kxf = 2tt/Lx]. This 
dependence has to do with the fraction of the energy in- 
jected by forcing that goes into the wavenumbers that 
contribute dominantly to the dynamo. In numerical ex- 
periments, increasing Lx does indeed lead to 7max ini- 
tially decreasing, but at large Lx it appears to stabilize at 
a constant value , a feature not captured by our the- 
ory. These simulations, however, had Rm = Re > 1, so 
the nonlinear mode coupling was nonnegligible there and 
the relation between the forcing and the velocity spectra 
less direct than in our calculation. In (astro)physically 



interesting systems, the approximation of linear shear is 
only valid at scales smaller than the characteristic system 
scale, which is rarely much larger than the turbulence 
scale, so, realistically, 27: / Lxkyf is never very small. 

Conclusion. We have presented what we think is the 
minimal analytically tractable model of the shear dy- 
namo. The simplicity of the model suggests that the 
effect is robust, while its rigorous validity in the assumed 
limits (Rm ^ Re ^ 1, Stc ^ 1) suggests that it is phys- 
ical. Much remains to be understood before it can be 
assessed whether the shear dynamo offers a panacea for 
(nonhelical) generation of large-scale magnetic fields in 
astrophysical systems. A further effort in this direction 
appears worthwhile in view of the great success recently 
enjoyed by shear-induced dynamos in astrophysically mo- 
tivated numerical experiments and a basic similarity of 
the field structure that they generate [13, 0, [HI . 
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